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Abstract 

Item banks, large sets of test items, can be used for the 
construction of achievement tests . Mathematical programming 
models have been proposed for the selection of items from an 
item bank for a test. These models make automated test 
construction possible. However, to find an optimal or even an 
approximate optimal solution to a test construction model can 
be time consuming. Here, it is shown how test construction 
models and in particular the Maximin Model are solvable by 
the program MPSX/370 V2. This program offers the user several 
implementations of the branch-and-bound method, which can be 
used £or solving test construction models. Several 
implementations are compared. The results show that tes^ 
construction models are solvable in a practical amount of 
time if the user applies the options offered by the program 
in an intelligent way. 

Keywords: Test Construction; Item Banking; Branch-And-Bound. 
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Implementations of the Branch-and-5ound method 
for Tast Construction Problems 

In this paper it is assumed that an item bank is available. 
The items of the bank are supposed to be calibrated under an 
item response model and some characteristics of the items 
such as format and content are supposed to be known. One of 
the applications of such an item bank is the selection of 
items for an in some sense optimal test, for instance, 
optimal with respect to its information function. 

An example of an item response model for binary scored 
responses is the 3-parameter model (Birnbaum, 1968) . In this 
model the probability, 2±{Q) , that an examinee with ability 6 
answers item i correctly is given by the formula 



Pi(G) = Ci + 



1-Ci 



1 + exp(-a i (6-b i )) 



where a^, b^ and c^ are the discrimination, difficulty, and 
guessing parameter of item i. The information function of 
item i is represented by 



ii<e) = 



a i 2 (l-c i ) 



(c i +exp(a i (6-b i ) )) (1+exp (~a L (Q-b L ) )) 2 



Another popular item response model is the Rasch model 
(Rascn, 1960) . In the Rasch model there is no guessing 
parameter (Cj_ = 0) and the item discrimination parameters are 
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equal. For a test with n : terns the test information function 
is found by addition of the item information functions: 

n 

1(6) = Z Ii(6). 
i=l 

The test information for an unbiased estimator of ability is 
the reciprocal of the (asymptotic) sampling variance of the 
estimator. Therefoxe, the higher the information function 
value, the better the test measures at a given ability level. 
This feature will be used in the construction of tests. A 
more detailed description of item response models and 
information functions is found in the psychometric literature 

(e.g. Lord, 1980; Hambleton & Swaminathan, 1985). 

The manual construction of nests from an item bank is 
almost practically impossible, because the number of possible 
tests is very large. Therefore, linear programming models 

(e.g. Theunissen, 1985; van der Linden & Boekkooi-Timminga, 
1989; Adema & van der Linden, 1989) are proposed which 
provide che test constructor with a computerized method for 
designing optimal tests. 

Linear programming (LP) is an optimization method 
applicable for the solution of problems in which the 
objective function and the constraints appear as linear 
functions of the decision variables (Rao, 1985) . In the LP 
models for test construction all cr almost all decision 
variables are of the 0-1 type. It is well known that in 
general models with integer variables are hard to solve. This 
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paper is addressed to the problem of solving test 
construction models, in particular the Maximin Model as 
proposed by van der Linden and Boekkooi-Timminga (1989) . 

In the following a short outline of the Maximin Model is 
given. Define the decision variables x^ as: 



x i = 



0 item i not in the test 



1 item : in the test, 



i ~* 1, 2 , • • * • / X, 



where I is the number of items in the item bank. The 
information function of the test, to be constructed is only 
considered at the ability levels 6 k , k = 1, . . . , K . The test 
constructor can specify the required precision by choosing 
the number and spacing of the ability levels b k . Let Ii(8 k ), 
i = 1, . «*, I be the information of item i at ability level 
0^. The relative amount of information required at ability 
level 0^ is specified by r k . The vector {r k } constitutes a 
target for the test information function. Let y be a decision 
variable such that (rtf, r K y) is a series of lower 

bounds to the information function of the test to be 
constructed. The idea of the Maximin Model is to maximize y, 
which implies that the lower bounds are maximized (maximin) . 
If n is the number of items to be selected for the test, then 
the Maximin Model can be written as follows: 



(1) Maximize y, 
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subject to 



I 



(2) 



2 Ii(e k )Xi - r k y £ 0, 



k = 1, 2 



/ • • • • / 



I 



(3) 



(4) 



Ax = b, 



(5) 



Xi € {0,1}, 



i a 1, 2 



/ • • • • / 



where x is (x^, X2/ . . . / xj) T , A is a mxl matrix, and b is a mxl 
vector. Constraints (2) considers the lower bounds to the 
test information function. Constraint (3) assures that the 
number of items in the test is equal to n. The constraints in 
(4) are added as a general provision to deal with practical 
constraints, for instance, on test composition, 
administration time, and the like; for examples, see Adema 
and van der Linden (1989). 

The purpose of this paper is to show how the Maximin 
Model caxi be solved by the computer program MPSX/370 V2. 
MPSX/370 V2 is an IBM licensed program for the handling of 
linear and mixed integer linear programming problems (LP 
problems with in-eger and continuous variables) . The program 
offers the usei a number of options for controlling the 
optimization process. The user is also provided with 
algoritmic tools which enable the experienced user to build 
his/her own heuristics and algorithms. A popular method for 
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solving mixed integer LP problems is the branch-and-bound 
(BAB) method. Several implementations of the F 3 method are 
available in MPSX/370 V2, which gives the user the 
opportunity to compare these implementations for test 
construction models . Several studies describing comparisons 
of implementations (strategies) in a more general perspective 
can be found in the mathematical programming literature (e.g. 
Breu & Burdet, 1974; Benichou et al. , 1971; Benichou, 
Gauthier, Hentges & Ribiere, 1977; Crowder, Johnson & 
Padberg, 1983; Ibaraki, 1987) . The importance of studies like 
these is the large amount of CPU-time that can saved by an 
effective implementation of the SAB method. In our case, fast 
solution methods are needed, for instance, to enable the test 
constructor to construct a test interactively. 

In the next section a general outline of the branch-and- 
bound method is given* More details about practical 
implementations of the BAB method are given further on in the 
paper. 



In this section the branch- and-bound method as described by 
Papadimitriou and Steiglitz (1982, pp. 433-438) is briefly 
reviewed for the case that the integer variables are of the 
0-1 type. Al_ the ideas to be presented carry over unchanged 
to the mixed 0-1 linear programming problem. Generally, a 0-1 
linear programming (0-1 LP) problem can be written as: 
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Problem 0: Max, z = c T x = c (x) , 
subject to 
Ax £ b, 

Xi e {0,1}, i = 1, , I, 



where c T = (c^, Cj) is the vector of cost coefficients. 

The branch-and-bound method starts wi^h the computation of 
the optimal solution x^ to the LP relaxation, i.e., Problem 0 
with 0 £ £ 1, i = 1, . I instead of e {0,1}, This 
computation is done by the simplex method, a well known 
method for sol/ing LP problems. The objective function value 
zq = c(x^) is an upper bound on thv* optimal objective 
function value z* - c(x*), where x* is the optimal solution 
to the 0-1 LP problem. If x^ is a 0-1 solution, then x* = x& 
and the 0-1 LP problem is solved. If x^ is not a 0-1 
solution then two subproblems are created. Suppose that, for 
instance, decision variable Xj in x^ is not equal to 0 or 1, 
Then the two subproblems are as follows: 



Problem 1: Max. z = c T x = c(x), 
subject to 
Ax £ b, 

Xj[ e {0,1}, i = 1, , I, 

x-i = 0, 
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Problem 2: Max. z = c T x = c(x), 
subject to 
Ax £ b, 
Xi e {0,1}, 
Xj = 1. 



X — 1 , a a a a , X, 



The variable Xj is called the branching variable. 

Both subproblems are solved by the simplex method- Then, 
one of the subproblems, say Problem 1, is chosen for 
branching. Several choice criteria are given in one of the 
forthcoming sections- Again a fractional decision variable, 
say x k , in x 1 is chosen and the Problems 3 and 4 are created 
by adding the constraint x^ - 0 and x^ = 1 to Problem 1. This 
process can be continued such thc*t a tree like in Figure 1 is 
created. 



Insert Figure 1 about here 



The root of the tree is formed by the LP relaxation and the 
children of a node are formed b} setting a fractional 
variable at its lower or upper bound (D or 1) . 

The branching process cannot ccnti - le indefinitely, 
because the number of variables is bounded- This implies that 
an optimal solution to the original 0-1 LP problem will 
always be found. Branching from a node is not needed in the 
following cases: 
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1) The LP solution corresponding to the node is a 0-1 
solution. Branching from this node is not needed, because 
adding constraints will give lers optimal solutions • If the 
objective function value of this 0-1 solution is better than 
che objective function value (z m ) of the best 0-1 solution 
obtained so far (incumbent) , then the solution becomes the 
incumbent and z m is now the objective function value of the 
new incumbent. If no 0-1 solution is known wh*;n the standard 
BAB method is started then the objective function value of 
the incumbent is initialized by z m = 

2) There is no feasible solution for the LP problem 
corresponding to the node, 

3) The objective function value, say z^ = cfxfy of the node 
is smaller than or equal to z m . This implies that a solution 
x that would be obtained as a descendant ui' x k would have 
objective function value c(x): 

C(X) £ Zfc <t z m . 

If there is no leaf left to branch on, the method stops and 
the ptimal solution to the 0-1 LP problem is the incumbent 
if one exists. 

The next two sections address a number of 
implementations of the BAB method. 
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Limiting the search 



The BAB method as given in the preceding section can be 
implemented in several ways. The following three ways of 
limiting the search are considered in this paper (see Adema, 
Boekkooi-Timminga & van der Linden, submitted) : 

LIM1: After the LP relaxation (Problem 0) is solved, 
variables are fixed using the reduced costs (see, e.g., 
Murtagh, 1981, p. 25) . Suppose d^ is the reduced cost 
corresponding to item i. The next two rules are used for 
fixing variables (H]_<1 is prespecif ied) : 

- Fix to 0, if in the relaxed solution = 0 
and zq - H]_zq < d^; 

- Fix x^ to 1, if in the ralaxed solution Xj_ = 1 
and zq - H]_zq < 

LIM2: Initialize z m by z m = &2 Z 0 for some prespecif ied H2<1 
instead of z m = -«>. This will decrease the number of 
branchings through the search tree. 

LIM3: Stop the search after the first 0-1 solution is found. 
This strategy is most useful in combination with LIM2, where 
H2 is close to 1, because this guarantees that the first 0-1 
solution found is a good solution. It is senseless to search 
for an exact solution, because the information functior 
values in the test construction models are stimatej and not 
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known exactly. 

The first two strategies speed up the BAB method 
considerably if H3 a/.d H2 are close to 1, because these 
strategies will reduc the number of variables and the search 
tree . An important characters 5tic of test construction 
problems is that in general the difference between zq a*»d z* 
is small, this characteristic makes it possible to choose H]_ 
and H2 close to 1. However, if H]_ and/or H2 are chosen to be 
close to 1, two kinds of problems can arise. Firstly, it is 
possible that too many varia' les are fixed at 0 and 1. This 
can imply that no feasible solution will be found. The second 
kind of problem occurs it r.o 0-1 solution exists with an 
objective function value between zq and ^2 Z 0- In b°th cases 
the branch-and-bound tree will be small and it does not take 
much time before it is clear that no feasible solution with 
objective function value higher than E2 Z 0 can ^ e found for 
the chosen values of H]_ and H2. The values of and/or H2 
can be adjusted and the solution procedure can be applied for 
the new values. A procedure for choosing H]_ and H2 will je 
given in the discussion section. 

Branching 

MPSX/370 V2 offers a variety of built-in strategies for the 
BAB method. The effectiveness oL a number of these strategies 
for test construction problems is regarded in the next 
section. In this section these strategies are described, but 
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first the notion of node estimation is incrocuc^ed. 

Node estimation 

A node estimation (IBM MPSX/370 V2 Program Reference 
Manual, 1988, p. 381) is a number that is attached to each 
node waiting for processing (waiting node) . J estimates the 
objective function value of the best 0-1 solution that can be 
expected at a descendent node and can be considered as a 
measure for the interest to continue the search by branching 
from the node. Here, two types of node estimations are 
considered. The first type of node estimation, say for node 
k, is found by computing 




I 

(6) E k = z k - 2 min [ PCL-^fi, PCUjMl-fj.) If 
i=l 



where 

E k = estimation for node k, 
z K = objective function value of node k, 
PCI*i = lower pseudo-cost of 0-1 variable x^, 
PCU^ - upper pseudo-cost of 0-1 variable x^, 

f^ = fractional part of the current value of 0-1 variable x^. 
Pseudo-costs were introduced by Benichou, Gauthier, Girodet, 
Hentges, Ribiere and Vincent (1971) and are well known in the 
literature of mathematical programming. They are used for 
predicting the detoriation of the objective function value 
when an integer variable with a fractional value is forced to 
take an integral value. During a search it is possible that 



17 



Branch-And-Bound 
13 

not all the pseudo-costs are known. Here, these missing 
pseudo-costs are assumed to be zero. For the computation of 
the pseudo-costs in MPSX/370 V2 the reader is referred to the 
IBM MPSX/370 V2 Program Reference Manual (1983, pp. 378-380). 

The second type of estimation for node Jc is found by 
computing: 

I 

(7) E k = - I mint f L , l-f L ) . 
i=l 

It should be noticed that the second type of estimation does 
not depend on the objective function value of node Jc. 

Branching Strategies 

The first branching strategy used in the forthcoming 
numerical experiments is: 

PURE: - Among the 0-1 variables with a fractional value the 
variable with the lowest subindex i is selected as 
branching variable. 
- Choice of the branching node: Let Jc be the last 
branching node. The next branching node is chosen 
according to the processing results of the two 
branches (k,n+l) and (k,n+2) (See Figure 2). 



Insert Figure 2 about nore 
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Three possible cases should be considered: 

1) Waiting nodes n+1 and n+2 satisfy the candidature 
conditions given in section "Branch-and-Bound for 
0-1 Linear Programming", i.e., there is no rule 
indicating that further branching from these nodes is 
senseless. In this case the node with the best 
estimation is chosen as branching node, where the 
first type of estimation (6) is used. 

2) Exactly one of the waiting nodes n+1 and n+2 
satisfies the candidature conditions. Then this node 
is chosen as branching node. 

3) Node n+1 and n+2 do not satisfy the candidature 
conditions. From among all the other candidate nodes 
the one with the best estimation (type 1) is chosen 
as the next branching node. 

The strategies PUREl through PURE4 differ from the PURE 
strategy only by the choice of the branching variable: 

PUREl: As branching variable the one with its value furthest 

from 0 and 1 is chosen. 
PURE2: As branching variable the one with its value closest 

to 0 or 1 is chosen. 
PURE3: As branching variable the one with its value closest 

to 0 is chosen. 

PURE4: As branching variable the one with its value closest 
to 1 is chosen. 
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In addition, the strategies INT through INT4 correspond 
to PURE through PURE4, the difference being that the choice 
of the branching node is based on the estimation of type 2. 
The last sec of strategies are USER through USER4 . These 
strategies correspond to PURE tnrough PURE4 . Now, however, 
the choice of the branching node is not based on node 
estimation but on the objective function value: The node with 
the highest value is chosen. 

The strategies PURE, INT1 , INT2 , INT3 , and INT4 are 
equal to the strategies with the same name in the IBM 
MPSX/370 V2 Program Reference Manual (1988) . In this manual 
the reader can also find a more detailed description of the 
strategies. 

Computational Experience 

The strategies described in the previous sections are 
compared in this section. The experiments were conducted with 
two simulated item banks containing 450 grammar items. The 
items in the banks fitted the Rasch {b^ _ N(0,1)) and 3- 
parameter model (a L _ U(0.5,1.5), bi _ Ni0,l), Ci = 0.2), 
respectively. 

In the Maximin Model as used in the experiments the 
ability levels were chosen to be $i = -1, 0 2 = 0, and 83 = 1. 
The relative information at these ability levels was set 
equal to 1 (r k =■ 1, k = 1, 2, 3). The Maximin Model was 
formulated as follows: 
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(8) Maximize y, 
subject to 

450 

(9) S Ii<6k> x i - Y * 0, k = 1, 2, 3, 
i=l 

450 

(10) L Xi - 20, 
i=l 

(11) x L € {0,1}, i a 1, , 450. 

Constraints (8) through (11) give the basic Maximin Model. As 
follows from constraint (10), the number of items in the test 
was required to be equal to 20. 

The experiments were conducted on an IBM9370. If an user 
wants to apply the options offered by MPSX/370 V2 in an 
advanced way, he/she has to write an ECL program for 
controlling the optimization process. The appendix shows the 
ECL program used in the experiments. ECL is a computer 
language based on PL\1. The CPU-times in the forthcoming 
tables are the times needed for the execution of the ECL 
program. 

In all the experiments strategies LIM1, LIM2, and LIM3 
were applied for limiting the search together with one of the 
branching strategies. The values of H2 in LIM2 were set to 
0.995 in all cases, implying that the difference between the 
objective function of the optimal 0-1 solution (z* = c(x*) ) 



9 

ERLC 



21 



Branch-And-Bound 
17 



and the 0-1 solution found was at most 0.5% of z*, because zq 
is an upper bOv^nd on z*. 

Table 1 and 2 display the CPU-times (in minutes) for the 
branch-and-bound part of the solution method and the 
objective function values for the computed 0-1 solutions. To 
be more specific the CPU-time is the time needed for the 
procedures MIXFLOW and SOLUTION in the ECL program (set 
Appendix) , because the CPU-times of those procedure depend on 
the branching strategy. In Table 1 the results for the bank 
with items fitting the Rasch model are shown. In this case 
two values for were used namely 0.995 and 0.9999. In Table 
2 the items fitted the 3-parameter model and was 0.995 and 
0.9975. Two values for were chosen to detect ne 

sensitivity of the results to the value. In the 3- 

parameter case the value can not be chosen as close to 1 
as in the Rasch case, because in the former case too many 
variables would be fixed, which would increase the chance of 
not finding a feasible 0-1 solution. 



Insert Table 1 and 2 here 



Looking at the two tables, it is seen that the differences 
between the CPU-times comparing the two values for are 
larger for the Rasch model than for the 3-parameter model. A 
possible explanation of this is the number of items fixed, 
which is much larger for the 3-parameter model for both 
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values of H^. 

The basic Maximin Model can be solved relatively easy, 
because the number of fractional values in the optimal 
solution for a relaxed (sub) problem is always smaller than or 
equal to the number of constraints . If a number of 
constraints is added to the model, a more interesting and 
realistic model is created. Therefore, the basic Maximin 
Model (8) -(11) was extended with the following constraints: 



450 

(12) S t^i £ 60, 
i=l 



150 

(13) I Xi = 7, 
i=l 



300 

(14) S Xi = 8, 
i=151 



450 

(15) S Xi - 5, 
i=301 



75 225 375 

(16) S Xi + 2 xi + I Xi = 12, 
1=1 i=151 i=301 



150 300 450 

(17) S Xi + S Xi + S Xi = 8. 
i=76 i=226 i=376 



Constraint (12) implied that tne administration time did not 
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exceed 60 minutes. The coefficients t± are estimates of the 
time a student from the population needs for answering item 
i. Here, they were drawn from the distribution U{1,6). The 
item banks were supposed to be partitioned in the following 
subsets with respect to content: 
Items 1-150: noun items; 
Items 151-300: verb items; 
Items 301-450: adjective items. 

Constraints (13) -(15) implied that the test contained 7 noun, 
8 verb, and 5 adjective items. Within the specified subsets 
the first 75 items were multiple choice items. e other 
items were of the matching type. According to constraints 
(16) and (17) the examinees had to answer 12 multiple choice 
and 8 matching items. In the experiments constraint (10) was 
left out of the LP model, because it was implied by the 
constraints (13) -(15) and (16) -(17). 

The same experiments as with the basic Maximin Model 
were conducted with the extended version. Table 3 and 4 show 
the CPU-times and the objective function values of the 
computed 0-1 solutions for the Rasch model and the 3- 
parameter model, respectively. 



Insert Table 3 and 4 here 



As in Table 1 and 2 the choice of the branching variable 
is important. In most of the experiments the best results in 
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terms of CPU-time were found when the branching variable was 
the variable furthest from an int*c>er (PURE1, INT1, and 
USER1) and closest to 1 (PURE4, INT4, and USER4) . Remarkably, 
the strategies INT through INT4 gave good results for all 
choices of branching variables when the \tems fitted the 3- 
parameter model but not when they fitted the Rasch model (see 
Table 1 and 3) . 

An interesting question is how hard is it to find an 
exact optimal 0-1 solution. Therefore, the basic and extended 
Maximin Model were also solved with the BAB method/ where 
LIM3 was not- applied, i .e . the algorithm was not stopped 
after the first 0-1 solution was found. LIM1 and LIM2 were 
applied with H]_ » H 2 , which guaranteed that an exact optimal 
0-1 solution was found if it had an objective function value 
higher than H 2 zq (see Discussion) . For the 3-parameter model 
Hi and H 2 were set equal to 0.995. For the Rasch model the 
gap between zq and z* is smaller, and therefore the values of 
Hi and H 2 were set equal to 0.9999. Only branching strategy 
PURE1 was applied. The CPU-times and optimal objective 
function values are displayed in Table 5. 



Insert Table 5 here 



The difference between z 0 and z* is indeed much smaller for 
the Rasch model. The reason is the larger similarity among 
the Rasch items, which makes it also harder to solve the 
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Maximin Model even if and H 2 are very close to 1. 

Discussion 

In this paper a comparison has been made among several 
search strategies for the branch-and-bound method. From the 
results it can be concluded that the strategies for limiting 
che search are very effective. Fixing a large number of 
variables by choosing H x very close to 1 and stopping the 
search after the first 0-1 solution has been found, reduces 
the CPU-time considerably. A good 0-1 solutior. is also 
guaranteed by choosing H 2 close to 1. Based on the 
experiments the following procedure for choosing the values 
Hi and H 2 is proposed; 

Ste » 1: If the items fit the Rasch model then = 0.9999 and 

H 2 - 0.995. if the items fit the 3-parameter model 

then = 0.995 and H 2 = 0.995. 
Step 2: Solve the Maximin Model with the BAB method. If a 

feasible 0-1 solution is found then go to Step 4 else 

go to 6tep 3. 

Step 3: If Hi > H 2 then H x := H 2 , else H 2 := H 2 *H 2 . Go to 

Step 2. 
Step 4: Stop. 

The starting values of ^ and H 2 are based on the results in 
Table 1 through 4. If Hj, m H 2 and no feasible 0-1 solution is 
found in Step 2, then no feasible solution exists with an 
objective function value higher than z m = H 2 z 0 (see, e.g., 
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Nemhauser & Wolsey, 1988, p. 389; Crowder, Padberg & Johnson, 
1983) , Hence, the adjustment scheme in Step 3 guarantees that 
the difference between the objective function value of the 
optimal 0-1 solution z* and the 0-1 solution found zp is at 
most 0.5%. If a smaller difference e is wanted, it can be 
realized by initializing H2 by H2 = 1 - e with e £ 0.005 and 
H l * H 2* Then / *>y adjusting the H]_ and H2 values according to 
Step 3 the algorithm in this paper is an e-approximation 
algorithm, i.e. (z* _ z F ) / z* £ e. 

Experiments with several branching strategies were 
conducted. The tables showed that especially the choice of 
the branching variable was important. Good choices for the 
branching variable were the variable furthest from an integer 
and the one closest to 1. The strategies INT1 through INT4 
are recommended in the IBM MPSX/370 V2 Program Reference 
Manual (1988) for models having a large number of 0-1 
variables as compared to the number of continuous variable 
and for which the user also tries to find a good solution but 
not necessarily the optimal. According to the tables these 
strategies were not much better, and in some cases even worse 
(see Table 1 and 3), if the Maximin Model is considered and 
the items fit the Rasch model. Table 5 shows that exact 
optimal 0-1 solutions can be found in a reasonable amount of 
CPU- time if LIM1 and LIM2 are used. For interactive test 
construction, however, searching for an optimal 0-1 solution 
can take too much time especially if the Rasch model is 
involved. 

The above conclusions are baseo -*n experiments with 
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models specified by the author. Nevertheless, they probably 
hold for Maximin Models in general, because the results 
showed the same tendencies for the branching strategies under 
varied conditions (with and without practical constraints, 
Rasch- and 3-parameter model) . The comparisons were also not 
influenced by the choice of computer implementation, computer 
characteristics etc, because all the experiments were 
conducted on the same computer and with the same program. 

In this paper it was shown how the computer program 
MPSX/370 V2 can be used for solving the Maximin Model. Other 
programs (e.g., HS/LP, Haverly Systems Inc.; SCICONIC, 
Sci conic Computer Services Ltd,; APEX4, Control Data 
Corporation) also provide the user with a number of options 
for controlling the search. Hence, the results can also be 
important for users of other optimization programs. The 
reader interested in computer codes for solving mixed integer 
linear programming problems is referred to Land and Powell 
(1979) and Powell (1985) . 



Appendix 



The program MPSX/370 V2 enables the user to control the 
search by an Extended Control Language (ECL) program, ECL is 
based on the computer programming language PL\1. The ECL 
program used for the experiments can be divided in three 
parts. In the first part the initialization takes place and 
the relaxed LP problem (Problem 0) is solved. Then, the 
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problem is reduced by fixing 0-1 variables at their bounds. 
Finally, the reduced problem is solved. 



/* NAME OF THE PROBLEM 



CONTROL: PROCEDURE OPTIONS (MAIN) ; 
/*****************************************************/ 

/* PART 1. INITIALIZATION AND OPTIMIZATION OF THE */ 

/* RELAXED PROBLEM */ 
/***************************************************** J 



%INCLUDE DPLINIT; 
DCL DUMMY $REAL; 
XBPNAME='BANK450' ; 
XDATA= ' BANK 4 50' ; 
XOLDNAME=' BANK45 0 ' ; 
CALL CONVERT; 
XMINMAX='MAX'; 
CALL SETUP; 
XOBJ=' OBJE' ; 
XRHS=' RIHA' ; 
CALL OPTIMIZE; 



CALL SOLUTION; 



/* OPTIMIZATION OF THE 
/* RELAXED PROBLEM 
/* PRINT THE RESULTS 



*/ 
V 
*/ 
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/A*********************************************** ***** f 

/* PART 2. FIXING OF 0-1 VARIABLES */ 
/*****************************************************/ 



XMXDROP=0.9975*XFUNCT; /* 


THE VALUE OF HI IS 


*/ 




/* 


SPECIFIED 


*/ 


nTTMMYr^YFrTMPT • 

U X — AC U IN X / 










, ' REVIFILE' , ' NOHEUR' , ' NOPRINT' ) ; 






1* 


THE VARIABLES TO BE 


*/ 




/* 


FIXED ARE FLAGGED 


*/ 


XPBNAME=' RBANK450' 


1* 


NAME OF THE REDUCED 


*/ 




1* 


PROBLEM 


*/ 


CALL REVISE ('FILE' 


, 'REVIFILE') ; 






/* 


THE PROBLEM IS REVISED 


*/ 




/* 


BY LEAVING OUT FIXED 


*/ 




/* 


VARIABLES 


*/ 



/****************************************************/ 

/* PART 3. OPTIMIZATION OF THE REDUCED PROBLEM */ 
/****************************************************/ 

CALL SETUP ( ' BOUND' , 9 BOUND' ) ; 



CALL OPTIMIZE; 

CALL MIXSTART; 
XMXMAXNO=l; 



XMXSTRAT='USER'; 



/* 
/* 

/* 
/* 
/* 
/* 
/* 



OPTIMIZATION OF THE */ 
RELAXED REDUCED PROBLEM*/ 

THE SEARCH STOPS AFTER */ 
XMXNAXNO 0-1 SOLUTIONS */ 
ARE OBTAINED */ 
THE SEARCH STRATEGY IS */ 
SPECIFIED BY THE USER */ 
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XMXSWT=' 0000004'; 


/* 


TO SPECIFY THE STRATEGY*/ 




/* 


7 SWITCHES ARE CHOSEN: */ 




/*' 


0000004=PURE1 */ 


XMXDROP=0 . 995*DUMMY; 


/* 


THE VALUE OF H2 IS */ 




/* 


SPECIFIED */ 


CALL MIXFLOW; 


/* 


SEARCH FOR 0-1 SOLUTION*/ 




/* 


USING THE BAB METHOD */ 


CALL SOLUTION; 


/* 


PRINT THE RESULTS */ 



END CONTROL; 

The information values I^Q^)/ k - 1/ 2, 3; i = 1, , 

were all multiplied by 1000 in the input file BANK450, 
because SETVAR did not work well with the non scaled values. 
The ECL program is given here to show how MPSX/370 V2 enables 
the user to control the optimization process and in 
particular how the optimization process can be controlled for 
the Maximin Model. For a description of the ECL language the 
reader is referred to the IBM MPSX/370 V2 Program Reference 
Manual (1988) . 
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Table 1 

Results for the branching cfrateaies applied to the basic 
Maximin Model where the items fit the Rasc h model. 



H 1 =0,995 a H]~0.9999 b 



Branching CPU-time Objective CPU-time Objective 
Strategy (mins . ) Function (mins „ ) Function 

Value Value 



PURE 


0.26 


3.9299 


0.09 


3.9299 


PURE1 


0.16 


3.9313 


0.07 


3.9313 


PURE2 


0.32 


3.9300 


0.08 


3.9300 


PURE3 


0.32 


3.9300 


0.07 


3.9300 


PURE4 


0.16 


3.9313 


0.06 


3.9313 


INT 


0.29 


3.9312 


0.16 


3.9312 


INT1 


0.24 


3.9287 


0.09 


3.9287 


INT2 


0.90 


3.9307 


0.15 


3.9307 


INT3 


0.89 


3.9307 


0.15 


3.9307 


INT4 


0.24 


3.9"37 


0.09 


3.9287 


USER 


0.25 


3.9299 


0.10 


3.9299 


USER1 


0.16 


3.9313 


0.06 


3.9313 


USER2 


0.33 


3.9300 


0.08 


3.9300 


USER3 


0.33 


3.9300 


0.08 


3.9300 


USER4 


0.17 


3.9313 


0.06 


3.9313 



a Number of fixed variables: 133. Elapsed CPU-time before 

MIXFLOW was called: 0.15 mins. 

b Number of fixed variables: 375. Elapsed CPU-time before 

MIXti,OW was called: 0.16 mins. 
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Table 2 

Results for the branching strate gies applied to the basic 
Maximin Model where the items fit the ^-par ameter model. 



H 1 =0,995 a H 1 =0.9975 b 



Branching CPU-time Objective CPU-time Objective 
Strategy (mins.) Function (mins.) Function 

Value Value 



PURE 


0.20 


4.i024 


0.16 


4.1024 


PURE1 


0.07 


4.1024 


0.04 


4.1024 


PURE 2 


0.24 


4.1043 


0.05 


4.1060 


PURE3 


0.25 


4.1043 


0.05 


4.1060 


PURE4 


0.07 


4.1024 


0.04 


4.1024 


INT 


0.08 


4.1033 


0.06 


4.1033 


INT1 


0.07 


4 .1043 


0.05 


4.1043 


INT2 


0.09 


4.1085 


0.06 


4.1085 


INT3 


0.09 


4.1085 


0.06 


4.1085 


INT4 


0.07 


4.1043 


0.05 


4.1043 


USER 


0.10 


4.1058 


0.09 


4.1058 


USER1 


0.06 


4.1024 


0.04 


4.1024 


USER2 


0.16 


4.0995 


0.05 


4.1060 


USER3 


0.16 


4.0995 


0.05 


4.1060 


USER4 


0.06 


4.1024 


0.04 


4.1024 



a Number of fixed variables: 393. Elapsed CPU-time before 

MIXFLOW was called: 0.14 mins. 

k Number of fixed variables: 419. Elapsed CPU-time before 

MIXFLOW was called: 0.15 mins. 
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Table 3 

Results for the branching strategies applied to the ^xtgnried 
Max.imin Moo^l wh ere the items fit the Rasch model. 





Hi- 


n qq ca 
u ♦ y yo a 




TT _ A 

H]_=U . 




Branching 


CPU-time 


0bj2ctive 




CPU-time 


Objective 


ocrauegy 


(ntins • ) 


Function 




(mins • ) 


Function 






value 






Value 


PURE 


0.23 


3.9308 




0.09 


3.9308 


PURE1 


0.25 


3.9295 




0.15 


3.9257 




1 7^ 






U ♦ Z 1 


o • yzu o 


rUKfciO 


u • DU 


O QO A Q 




fi 1*7 

U • 1 / 


"3 QIAO 






Q QQ A Q 




n a q 

u ♦ uy 


*3 QOrtO 

o • y oUo 


INT 


0*43 


3.9306 




0.19 


3.9306 


INT1 


0.S9 


3.9296 




0.19 


3.9296 


INT2 


2*59 


3.9245 




0.17 


3.9299 


INT3 


1.91 


3.9310 




0.14 


3.9310 


INT4 


0.38 


3.9212 




0.14 


3.9212 


TTQPP 


0.22 


3.9308 




n no 
U . Uo 


O Q Q A Q 


USER1 


0.21 


3.9295 




0.14 


3.9257 


USER2 


1.35 


3.9308 




0.21 


3.9206 


USER3 


0.60 


3.9308 




0.17 


3.9308 


USER4 


0.22 


3.9308 




0.09 


3.9308 


a Number 


of fixed 


variables: 


133 


. Elapsed 


CPU-time 


MIXFLOW 


was called: 0.19 mins 








k Number 


of fixed 


variables: 


377 


. Elapsed 


CPU-time 


MIXFLOW 


was called: 0.20 mins 
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Table 4 

Results for the branching strategie s applied to t.ho extended 
Maximin Model where the jfcg mS fit the ^paramfiter model , 



H!=0.995 a H 1 =0.9975 b 



Branching CPU-time Objective CPU-time Objective 
Strategy (mins.) Function tmins.) Function 

Value Value 



PURE 


0.43 


4.0531 


0.17 


4.0583 


PURE1 


0.13 


4.0583 


0.20 


4.0607 


PURE2 


0.84 


4.0580 


0.37 


4.0607 


PURE3 


0.36 


4.0565 


0.48 


4.0580 


PURE4 


1.20 


4.0580 


0.14 


4.0607 


INT 


0.12 


4.0580 


0.09 


4.0580 


INT1 


0.09 


4.0589 


0.09 


4.0580 


INT2 


0.41 


4.0575 


0.18 


4.0575 


INT3 


0.41 


4.0575 


0.17 


4.0575 


INT4 


0.12 


4.0580 


0.08 


4.0580 


USER 


0.27 


4.0583 


0.17 


4.0583 


USER1 


0.25 


4.0583 


0.17 


4.0583 


USER2 


0.56 


4.0607 


0.78 


4.0589 


USER3 


0.30 


4.0565 


0.39 


4.0566 


USER4 


0.47 


4.0607 


0.26 


4.0583 



a Number of fixed items: 390. Elapsed CPU-time before MIXFLOW 

was called: 0.17 mins. 

b Number of fixed items: 414. Elapsed CPU-time before MIXFLOW 

was called: 0.17 mins. 
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Table 5 

Exact optimal 0-1 solutio ns by branching strategy PHRE1 

Objective 
Function Value 



IRT 


Extended 


Relaxed 


0 


-1 


CPU-time 


Model 


Maximin 








(mins. ) 


Rasch 


no 


3.9318 


3 


.9317 


3.24 




yes 


3,9314 


3 


.9313 


4.23 


3-param. 


no 


4,1188 


4 


.1105 


0.82 




yes 


4.0731 


4 


.0607 


2.64 
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Figure Captions 



Figure 1 . Branch-and-bound tree. 



Figure 2 . Node k and its children. 
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